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ABSTRACT 

We present Wide Field Spectrograph (WiFeS) integral field spectroscopy and HST EOS spectroscopy 
for the LINER galaxy NGC 1052. We infer the presence of a turbulent accretion flow forming a small- 
scale accretion disk. We find a large-scale outflow and ionisation cone along the minor axis of the 
galaxy. Part of this outflow region is photoionised by the AGN, and shares properties with the ENLR 
of Seyfert galaxies, but the inner (i? < 1.0 arcsec) accretion disk and the region around the radio 
jet appear shock excited. The emission line properties can be modelled by a “double shock” model 
in which the accretion flow first passes through an accretion shock in the presence of a hard X-ray 
radiation, and the accretion disk is then processed through a cocoon shock driven by the overpressure 
of the radio jets. This model explains the observation of two distinct densities 10^ and ~ 10® cm“®), 
and provides a good fit to the observed emission line spectrum. We derive estimates for the velocities 
of the two shock components and their mixing fractions, the black hole mass, the accretion rate needed 
to sustain the LINER emission and derive an estimate for the jet power. Our emission line model is 
remarkably robust against variation of input parameters, and so offers a generic explanation for the 
excitation of LINER galaxies, including those of spiral type such as NGC 3031 (M81). 

Subject headings: physical data and processes: black hole physics, shock waves - galaxies: individual 
(NGC1052, M81) -galaxies: jets - galaxies: nuclei 


1. INTRODUCTION 

Since the low-ionisation nuclear emission line region 
(LINER) class o f galaxies was originally formalised by 
Heckman (19801, the nature of their exci t ation has re- 
mained a cause of speculation. Heckman (1980) distin¬ 
guished these galaxies from Seyfert galaxies and from 
galaxies with H Il-region-like emission lines with the 
following criteria: [O II] A3737,9> [O III] A5 007, and 
[O I] A6300> [O III] A500 7. The su r veys b y |0’ Con- 
nell &: Dressel (1978) and Heckman ( 1980| confirmed' 
what had previously been suspected from a few coin¬ 
cidences: a strong correlation between bright LINER 
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emission and powerful compact nuclear radio sources. 
In most of these bright LINERs, virtually all of the 
radio emission comes from fla t-spectrum compact self - 
absorbed synchrotron sources (jCondon & Dressel||1978|). 
Strong LINER emission is also sometimes seen in pow¬ 
erful radio galaxies with exten ded steep-spectrum syn¬ 
chrotron em ission, such as M87 (Ford et al.|l994 Harms 
et al.|[T9M |. 


Low level LINER emission has now been detected 


in a l arge fraction of elliptical galaxies ([Phillips et al. 
19861 Veron-Cetty fc Veron||I986| |Ho||I996[ |Ho et ^ 


1997a r This emission may have a variety of origins 


nonetheless the bright LINERs associated with power¬ 
ful radio sources appear to be true AGN. Both their 
radio emission, and the presence of b road Ho emission 
(Ho et al. 1997b Younes et al. 2012) suggests a link to 
radio-powertul beytert galaxies and quasars. A few LIN¬ 
ERs have been observed to have strong, compact and 
spectrally-hard X-ray emission implying t hat the source 


is an A GN (Serlemitsos et al 1996). Indeed, Kewley et al. 
(2006) found that the nuclear LINERs in red emission¬ 
line SDSS galaxies lie along the same host-galaxy Ed¬ 
dington ratio relations as Seyferts, which also implies 
an AGN origin for these LINERS. Energy considerations 
suggest that the power source is the release of gravita¬ 
tional energy by matter dropping into a supermassive 
black hole. K inematic measurements, beginning w ith the 
LINER M87 ( Ford et al.||1994 Harms et al.||1994 ), show 
that the velocity widths can be understood as primarily 
due to orbital motion around a massive black hole. How¬ 
ever, even the “narrow line” components in LINERs have 
FWHM of hundreds of km s“^ when observed through 
HST apertures. This is broader than expected from the 
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variation of rotational velocity across a small aperture 
placed off-nucleus, and implies the presence of turbu¬ 
lence, non-rotational flows, or high velocity shocks. 

The optical emission line ratios of LINERs were origi 


nally modelled as being due to shock-heated gas (Bald- 


win et al.l 1981). It was soon discovered, however, that 


very similar liiie ratios could also be produced by gas 
that was photoionized by either a power-law or a thermal 
Bremstrahhlung continuum with a low ionization param- 
eter; i.e., a low ratio of i ncident ionizing photons per ion 
dFerlan d &: Netzer|1983D. To add to the possible physical 
models, TlHevicii &: JVlelmck ( 1985|) pointed out that a 
similar hard U V continuum could be produced by very 
hot luminous Wolf-Rayet stars in a nuclear starburst. 

Strictly speaking, the term LINER should be confined 
to describing the nuclear re gion of a galaxy - according 
to the original definition by Heckman (19801. However, 
many host galaxies show extended emission on kiloparsec 
scales which maintains a LINER-like spectrum. Many of 
the objects with LINER-like spectra found in surveys 
are of this nature, and it may well be that this extended 
emission is excited in a different manner from the nu¬ 
clear emission. Such regions were distinguished from the 
nuclear emission li ne regions as “ex t ended LINERs” and 
were analysed by Yan & Blanton (2012). They found 
that the extended emission is inconsistent with ionisation 
from a central object, and that rather, extended sources 
such as post-AGB stars are required to p rovide the ion¬ 
isation ; an idea originally proposed by [Binette et al.| 
(1994), and recently developed in t he light of CALIFA 


integral held spectroscopic data by Singh et al. (2013). 
These authors conclude that extended LliNER emission 
of itself cannot be use to infer the presence of an AGN. 

From the velocity dispersions and the small width dif- 
ferences between different lines, Yan & Blanton (2012) 
concluded that the line width in extended LliNERs is 
dominated by bulk motion and turbulence, and that two 
ISM components are needed to produce line width dif¬ 
ferences between different ions. This is perhaps difficult 
to reconcile wit h ionisation by the p ost-AGB component 
alone. However, Farage et al. (2010) showed that shocks 
in a multi-phase medium during a merger event can pro¬ 
duce different line widths i n different ions. Fur thermore, 
in hrst-ranked Ellipticals (Farage et al. 2012), the ob¬ 
served line ratios show (apparently) more shock domi¬ 
nance closer to the nucleus, again showing how difficult 
it is to interpret the extended LINER emission using a 
single mechanism. 

Significant progress in LINER emission line diagnos¬ 
tics can be made by extending the flux predictions for 
power law, bremsstrahl ung, and shock mod els into the 
UV. In the case of M87, Dopita et al. (1997) have shown 
that several UV lines are much stronger m shock mod¬ 
els than in the other types of model. Unfortunately, the 
UV spectral data is av ailable for relativel y few LINERs. 
Such data as do exist (Maoz et al. 1998) seem to show 
that the LINER class is not homogenous in its UV prop¬ 
erties, with some objects exhibiting strong emission lines, 
and others having a UV spectrum that is simply consis¬ 
tent with an old stellar population. 

NGG 1052 is an E4 galaxy, is the prototypical LINER 
galaxy, and being nearby (z= 0.0045) can be studied in 
detail. It is a radio galaxy having a well-c ollimated non- 


thermal radio jet on both arcsec scales (Wrobel 1984 



R.A. Offset (mill) arc sec.) 



R.A. Offset (arc sec.) 


Fig. 1.— Upper panel, from|Claussen et al.|||1998|l shows the ra¬ 
dio jet and the H20 masers in iMldU iUt>2. Une mas corresponds to 
a distance of about 18000 AU. The masers are probably associated 
with a dense circum-nuclear disk, and are formed by amplifica¬ 
tion of the background radio continuum emission. The lower panel 
shows the HST image of the (galaxy subtracted) Ha jet - perfectly 
aligned with the inner radio structure, but about 500 times longer. 
This structure apparently represents the interaction of an outflow 
with the circum-nuclear ISM. The white regions are regions of high 
dust extinction and may delineate the shock-compressed cocoon of 
the radio jet. At larger distances (~ 4 arcsec) from the nucleus 
there are a number of bright knots, which are par t of a filamentary 
complex best seen in fig. 2 of|Pogge et al.| ||2000|l. Figure courtesy 
of Dr. Andrew Wilson, deceased. 


Cooper et al. 


et al.||1998D. T 


2007) and milli-arcsec scales (Claussen 


Sion with 1 3 

eTaLjlMS ). 


re radio jet s hows sub-relativistic expan- 
0.24 — 0.38 (Kadler et al. 2004b Lister] 


However, it also has a dense nuclear accretion disk, 
demonstrated by the fact that it is one of only two el¬ 
liptical galaxies containing very luminous water masers. 
The H 2 O maser emission lines are relatively broad 
(90 km s“^) and smooth, unli ke the narrow mase r spikes 
seen in other ac tive galaxies ( Br aatz e t al.||1996). VLBI 
observations by |Claussen et m.| (|l998|) (see hg. Tj) show 
that the masers lie along, rather than perpendicular to, 
the jet, and only on the west south-west side of the nu¬ 
cleus. The sc ale of the maser emiss ion is about 400 /ras 
(^7000 AU). iClaussen et al.|(|I998|) argue that although 
the jet is energetically capable of powering the observed 
maser emission by driving slow, non-dissociative shocks 
into circum-nuclear, dense molecular clouds, it is more 
likely, given their location and velocity structure, that 
the masers represent the direct amplification of the radio 
continuum emission of the jet by the foreground molec¬ 
ular clouds in an inclined nuclear accretion disk. 

At X-ray wavelengths, the nucleus appears as a hard, 
heavily obscured (Vh « 10^^ cm“^), flat-spectrum (a ~ 
0.1) point source with an (unabsorbed) luminosity of 
5 X lO^^erg s“^ (Guainazzi & Antonelli 1999 Kadler 
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et al.|2004b )■ With CHANDRA observations, a more ex- 
tended region of thermal X-ray emission {kT ^ 0.5 keV) 
is seen occupying a double wedge-shaped region extend¬ 
ing 10 ar csec. east and west of the nucleus (Kadler 
et al.||2004b|). Alt hough it s exten t is simila r to the radio 
lobes o bserved by |Wrobel| (|1984|), in detail, |Kadler et al. 

(2004b I finds an anti-correlatidn between tne X-ray and 
radio surface brightness, suggesting that the X-ray gas 
may be associated with the working surface of the radio 
lobes. 

The optical emission line gas (see fig. shows a jet 
component extended on arcsec scales, perfectly aligned 
with the innermost radio jet, possibly ionised by UV syn¬ 
chrotron emission beamed from the inner jet. Elsewhere 
there are Ha emitting blobs and filaments filling the zone 
where the thermal X-rays are observed. However, the 
most prominent feature as w ell as a bright nuc lear re¬ 
gion 100-200 pc in diameter (j Pogge et al.|[2000l). This 
has been studied in detail by |Sugai et al. Ip2005|), and 
probably represents a resolved accretion disk. I’he line 
flux coming from the much larger extended region is sim¬ 
ilar in total luminosity to that of this inner disk. 

In some sense, NGC 1052 could be classified as a 
“hidden broad-line LINER” following the detection of 
a prominent broad Ha com ponent with EW HM ^ 
5000 km s“^ in polarised light (Barth et al.|l999 ). Broad 
wings are also seen in un-polarised light. 'I'his data es¬ 
tablishes the existence of a dense inner accretion disk 
with a hot scattering region around it. 

The interpretation of the emission line spectrum of 
NGC 1052 (like other LINERs) remains controversial. 
The idea that the NGC 1052 emission is shock excited 


was originally suggested by Koski & Osterbrock (1976), 
and a detailed radiat ive shock model was presented by 
(Eosbury et al. 1978). However, this model fell out of 
favour following the apparent success of photoionisation 
model s. However, from infrared data, |Sugai fc Malkan] 
(2000) found evidence that shock excitation might well 


provide a substantial fraction of the total emission in 
NGC 1052. In contradiction to this explanation, a n ex¬ 
amination of the HST data by Gabel et al. (2000|) (re¬ 
analysed here) suggested that the the emission-line huxes 
can be simulated with a simple photoionization model us¬ 
ing a central power-law continuum source with a spectral 
index of a = —1.2. Two components with radically dif¬ 
ferent densities and ionisation parameters were required 
in this model, leaving the physical origin of these two 
components something of a mystery. 

A vital clue to the resolution of this issue was provided 
by the sub- arcsecond reso lution integral field observa¬ 
tions of Sugai et al. (2005) who used the Kyoto Tridi- 
mensional Spectrograph 11 mounted on the Subaru Tele¬ 
scope. They found evidence for a high-velocity bipolar 
outflow, a low-velocity rotating disk and a spatially unre¬ 
solved broad line region visible in H/3. The high velocity, 
outflowing [O III] emission ridges closely correlate with 
th e radio structures ma pped wit h MERLIN at 1.4G Hz 
by |Kadler et al.| (|2004a|) . Whilst |Sugai et al.| (|2005|) in- 
terpret these observations in terms of a directed outflow, 
we argue here that it is more likely that this emission 
is caused b y a cocoon shock as advocated by |Bicknell| 
et al. (1997) in the context of gigahertz-peaked spectrum 
(GHSj yrd compact steep-spe ctrum (CSO) radio sources 
and by Bicknell et al. (1998) in the context of Seyfert 


galaxies. 

In this paper we argue that, in NGC 1052, such cocoon 
shocks are propagating into the (already shocked) accre¬ 
tion disk, giving rise to the dense emission component re¬ 
quired by the HST observations. In addition, we present 
integral field spectroscopy of the central 25 x 38arcsec 
region of NGC 1052 which establishes the existence of 
a higher excitation “ionisation cone” and a large-scale 
bipolar mass outflow presumably energised by the radio 
jet. In addition we will show the existence of a turbulent 
accretion flow along the major axis of the galaxy, and 
in the plane of rotation of the stars. In order to under¬ 
stand the excitation of the i nner, bright Ha reg ion, we 


re-examine the HST data of Gabel et al. (|2000 ), in the 


light of the observations of Sugai et al. ) and con¬ 


struct a new accretion shock plus cocoon shock model 
which naturally explains the existence of the two density 
components referred to above. 

The paper is organised as follows. In Section we 
present the data, and the extraction of the reduced quan¬ 
tities. In Section!^ we introduce the physical basis of 
our multiple shocK model to explain the spectrum of 
NGC 1052. In Section]^ we present our accretion -|- co¬ 
coon shock model in the presence of a hard X-ray spec¬ 
trum which is very successful in reproducing the observed 
LINER spectrum of NGC 1052 over a wide range of as¬ 
sumed shock velocities in either component. In Section]^ 
we use the observations to estimate the black hole mass, 
the mechanical energy luminosity of the radio jet, and the 
mass accretion rate into the accretion disk. We assume 
the redshift-independent distance to NGC 1052 given 
by the NASA Extragalactic Database (NED); 19.9 Mpc 
which implies a spatial scale of 97 pc arcsec”^. 

2. OBSERVATIONS AND RESULTS 
2.1. The WiFeS integral field observations 

NGC 1052 was observed on 2-3 No vember 2013 using 
the W ide Eield Spectrograph (WiEeS; Dopita et al.|2007 


2010) at the ANU 2.3m telescope at Siding Spring Ob- 


servatory. WiFeS is an optical integral field spectrograph 
providing a field of view of 25 arcsec x 38 arcsec via a 
total of twenty five 38 arcsec x 1 arcsec slitlets. The in¬ 
strument is a double-beam spectrograph, with arms for 
both the blue and for the red part of the spectrum. Thus, 
a wide wavelength range in the optical is covered simul¬ 
taneously, and with adequate spectral overlap between 
the red and blue spectra. 

The observations for NGC 1052 consisted of 3 x 1000s 
exposures at RA = 02:41:04.9 Dec= -08:15:21.1 (J2000) 
and 3 x 1000s exposures on the nearby sky at the position 
RA = 02:41:10.5 Dec= -08:11:41.4 (J2000). All observa¬ 
tions were obtained using the Rs = 3000 grating in the 
blue (733000 grating) and the Rs = 7000 grating in the 
red arm (777000 grating). The position of the WiFeS data 
cube is shown on an image of NGC 1052 in Figurej^with 
the physical scale at the assumed distance of 19.5 kpc. 
Absolute photometric calibration was made using the 
STIS spectrophotometric standard stars HD 009051 and 
HD 031128 In addition two B-type telluric standards 
were observe, HIP 108975 and HIP 18926. During the 

Available at : 

WWW.mso.anu.edu.au/ bessell/FTP/Bohlin2013/G012813.html 
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observation, the seeing averaged 1.3-1.5 arcsec, reason¬ 
ably well matched to the 1.0 arcsec pixels of the spectro¬ 
graph. 
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Fig. 3.— The 38 x 25 arcsec WiFeS image of NGC 1052 in [O III] 
A5007 (blue), Ho (green) and [N II] A6584 (red). The stretch has 
been chosen so as to enhance the fainter features, and the image has 
been smoothed with a rectangular box-car function so as to reduce 
pixelation. Note the E-W extension of the [O III] - bright c entral 
region al ong the direction of the arc sec. scale radio jets ( fWr^ 
|bel|11984|l, and the prominent X-shaped central ionisation cones m 
[u 111]. Ho; is enhanced in a two-sided funnel in the direction of 
the minor axis of the galaxy (along the rotation axis of the stars). 
[N II] is enhanced along the major axis in the region of extended 
thermal X-ray emission. 


Fig. 2.— The SDSS g,r,i image of NGC 1052 with the WiFeS 
pointing position superimposed. We also show the scale at the 
assumed distance of 19.9 Mpc. 


2.1.1. WiFeS Data reduction 

The data were reduced using the PyWiFeS pip eline 
written for the instrument (Childress et al.|2014|. In 
brief, this produces a data cube which has been wave¬ 
length calibrated, sensitivity corrected (including telluric 
corrections), photometrically calibrated, and from which 
the cosmic ray events have been removed. 

For the resultant data cube, we have used the IFS 
toolkit LZIFU (Ho et al. 2014 in prep.) to extract 
gas and stellar kinematics from the WiFeS data. LZ- 
IFU uses the the penalized p ixel-fitting routine (ppxf 
Cappellari & Emsellem 2004) to perform simple stel- 
lar population (SSF) synthesis fitting to model the con¬ 
tinuum, and fits the emission lines as Gaussians us- 
ing the Levenber g-Marquardt least-squares technique 
(Markwardt 2009). For generating stellar velocity and 
velocity dispe rsion maps we employ the theor etical SSP 
libraries from Gonzalez Delgado et al. (2005) assuming 
the Padova isochrones, h'or the gas velocity and velocity 
dispersion maps we simultaneously fit within each spaxel 
one-component Gaussians to 11 strong optical emission 
lines and we constrain all the lines to have the same ve¬ 
locity and velocity dispersion. The 11 lines of the fit are 
[O 11]AA3726,29, H/3, [O 11I]AA4959,5007, [O 1]A6300, 
[N 11]AA6548,83, Ha, and [S II]AA6716,31. We fix the 
ratios [O III]AA4959/5007 and [N ll]AA6548/6584 to their 
theoretical values g iven by quantum mechanics (|Dopita 
' fc Sutherland||2003 ). 


the degree of excitation is in no sense uniform. Indeed, 
in a restricted X - shaped zone extending nearly 10 arc¬ 
sec from the nucleus, we find [O III] lines to be strongly 
enhanced, approaching the values seen in Seyfert narrow 
line regions; see figure We regard this as probably 
delineating an ionisation cone photo ionised by the cen¬ 
tral engine. In the nuclear region an enhanced region 
of [O III] emission i s seen to be a ssociated with the arc 
sec. scale radio jet ( Wrobel||1984), while Ha is enhanced 
along the axis of rotation of the stars (the minor axis). 
The [N II] emission is relatively strongest along the major 
axis of the galaxy, within the X - shaped [O III] ionisa¬ 
tion cone in which the [O III]/H/3 ratio rises to 3-5. The 
spatial extent and shape of this region is very similar to 
the region within which th e thermal X-rays are enhanced 
(c./|Kadler et al.| (|2004bD). 

The dynamical picture - shown in Figures]^- is equally 
interesting. The stars rotate smoothly around the pho¬ 
tometric minor axis, and are characterised by a velocity 
dispersion of » 200 km s~^, not unusual for an ellip¬ 
tical galaxy (|Ge rhard et al. 2001 Thomas et al. 2007 


Gappellari et al.||2013p. There is a sharp cusp in the ve- 


2.1.2. Results from the Data Cube 

The WiFeS data cube furnishes a number of novel and 
interesting results. The emission line maps reveal that 


locity dispersion close to the nucleus, presumably within 
the zone of influence of the central black hole. The dy¬ 
namics of the gas is completely different to the stellar 
dynamics. The gas velocity component analysis reveals 
two bubbles of gas expanding along the minor axis of the 
galaxy, separated by a region of very high velocity dis¬ 
persion. These bubbles are about 15arcsec in diameter 
(^ 1.5 kpc), and are well-aligned with the Ha emission 
“funnel” seen in Figure The bubble on the NE side 
is approaching, so this is located on the near side of the 
galaxy. This is consistent with the fact that the H 2 O 
maser emission associated with amplification of jet radio 
continuum by the inner accretion disk is on the opposite 
side of the nucleus, also sugge sting that the ENE j et is 
directed somewhat towards us (Glaussen et al.|19'M) (see 
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Fig. 4.— The dynamical structure if NGC 1052 in Ha and in the stellar component. The Ha velocity field (panel a.) shows two bubbles 
of gas expanding along the minor axis of the galaxy, separated by a region of very high velocity dispersion seen in panel b. The stellar 
dynamics are dynamically quiescent. The stars rotate slowly around the photometric minor axis (panel c.), and with a velocity dispersion 
of 5^ 200km s“^, typical of the values observed in elliptical galaxies (panel d.). 


fig- 0- The velocity of approach in the NE bubble, and 
the recession velocity in the SW bubble is greatest along 
the axis, but decreases strongly towards the boundary 
of the bubbles. In addition, the internal velocity field 
is irregular and quite unlike a rotation curve. Further¬ 
more, both bubbles show an internal velocity dispersion 
of about 120 km s“^, much greater than the instrumental 
resolution in the red (~ 45 km s“^. All these dynamical 
properties are consistent with the buoyant expansion of 
turbulent bubbles filled with hot or non-thermal plasma 
in the ISM of the galaxy. 

The region between the two expanding Ha-emitting 
bubbles displays a very high velocity dispersion, up to 
300 km s“^. In both extent and morphology it coincides 
with th e region within which the thermal X-rays are en¬ 
hanced (Kadler et al.|2004b ) and the region in which the 
[N II] emission is unusually strong with respect to Ha, 
referred to above. The juxtaposition of high velocity dis¬ 
persion, enhanced [N II]/Ha ratios, and the presence of 
thermal X-ray emission all suggest the presence of str ong 
radiat i ve sho cks in a highly turbulent medium c.f |Ho| 


et al. (2014). Its inner region also correspcmds to the 
zone of enhanced obscuration seen in Figure [l] 


2.2. The FOS Observations 

The spectral data used in this study were obtained un¬ 
der project GO-6532 (investigators: Dopita & Koratkar) 
using the Faint Object Spectrograph (FOS) aboard the 
Hubble Space Telescope (HST /FOS) on 1997 January 13 


(post-COSTAR). All observations were made using the 
0.86 arc sec. pair square aperture and the wavelength 
coverage was 1200 - 6800A. The optical spectra (G400H 
and G570H grating) and the NUV (G270H grating) data 
were obtained at a resolution of A/AA ~ 1300 (or an 
approximate resolution 2.1-4.4A). Because of the faint¬ 
ness of the lines, the UV spectrum was obtained with 
the G160L grating giving a resolution of A/AA ~ 250 (or 
an approximate resolution of 6.4A). For NGC 1052, the 
aperture was well-matched in size to the central bright 
region visible in the Ha image of Fig. and so provides 
the integrated spectrum of the true nuclear LINER in 
this galaxy. 

The G160L data were taken in three exposures cover¬ 
ing two and a half orbits. These were followed by the 
G270H exposure in the next orbit and the G400H and 
G570H exposures in the final orbit. Only one guide star 
was acquired at the outset, so the roll of the telescope 
was controlled by gyroscope, which generally results in 
a slow drift of the target across the detector over the 
course of several orbits. Archival data for the two full- 
orbit G160L exposures were checked for systematic dif¬ 
ferences in the bright emission lines that would indicate 
substantial drift. Figure shows that such differences 
are small. 


2.2.1. Extracting the Emission line spectrum 
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Fig. 5.— Archival spectra from two orbit-long G160L exposures 
(solid and dashed line) and their difference (dot-dashed line), com¬ 
pared for evidence of drift of the target. We may conclude that 
any such drift is negligible. 


Figureshows the FOS spectrum, along with the con¬ 
tinuum and line model fits used to extract the line fluxes. 
We use a modified version of lzifu to fit the spectrum. 
For the continuum fitting, w e adopt the MILES stellar li¬ 
braries (Vazdekis et al. 20101 at wavelengths greater than 
A « 3500A. At wavelengths smaller than A « 3500A, 
we fit the continuum with linear combinations of 10 or¬ 
ders of Legendre polynomial. Channels around emission 
lines are masked out while fitting the continuum. For the 
emission line fits, we split the emission lines in 14 groups 
based on their wavelengths and then perform separate 
fits to each group. In each fit, we constrain all the lines 
to have the same velocity and velocity dispersion. All the 
emission lines are described by single component Gaus- 
sians, except for the Ha line where a second Gaussian 
component is required to properly describe the line pro¬ 
file. Where possible, the flux ratios of the emission lines 
from the same multiplet of a given ion are constrained 
based on the relative probabilities of their atomic tran¬ 
sitions. 


2.2.2. The FOS Emission Line spectrum 

The de-reddened extracted emission line intensities are 
listed in Table 1. The line fluxes have been de-reddened 
by the the oretical function for a phy sical grain model 
derived by Fischera & Dopita (2005) and which pro¬ 
vides a very close fit to the Galactic attenuati on law with 
R = 3.2. In common with Gabel et al. ([2000) we find an 
extinction which is unusually high tor the LINER popu¬ 
lation; Ay = 1.05. This corresponds to E(B-V) = 0.33, 
somewhat lower than the val ue of E(B-V) = 0.44 esti¬ 
mated by [Gabel et al.| ([2000| . The presence of such a 
large amount of circuin-nuclear dust suggests that the 
nuclear region also contains a large amount of interstel¬ 
lar gas, which we will interpret as due to a dusty ac¬ 
cretion flow onto the bright nuclear disk. Inde ed, this 
dusty flow ca n be seen directly in Figure 4 of [Kadler[ 
et al. (2004a) as a complex of accretion streams within 
the central ~ 2 arcsec 

The line fluxes of Table 1 cannot be fit by a simple 


one-zone ph otoionisation model. As realised by Gabel 
et al. (2000) the line ratios are inconsistent with such 


a model. While the [S II] A6717,31 line ratio and the 
presence of a strong [O II] A3726,9 doublet suggest a 
density of order 300 < Ue < 1000 cm“^, ratios such 
as [S II] A4068,76/6717,31, [N II] A5755/6584, [O II] 
A2470/3726,9 and [O III] A4363/5007 are all ind icative 
of densities in the range 10® 10^ cm“® 


Gabel 


et al. ([2000[) model these two components as physically 


independent entities photoionised by the same power- 
law spectrum, and derive appreciably different ionisa¬ 
tion parameters for the two components. Surprisingly 
the denser component has the higher ionisation parame¬ 
ter, so if the model is correct, it would have to be emitted 
much closer to the central engine. If so, why then does 
the dense component not show the broader line widths 
expected with the Keplerian motion around the central 
engine? In what follows we will present a physically mo¬ 
tivated accretion model with photoionisation which ap¬ 
pears to provide a natural answer to this question. 


3. MULTIPLE SHOCKS IN NGC 1052 
3.1. The Large Scale Flows 

The WiFeS observations of NGC 1052 clearly reveal 
the presence of a high velocity dispersion gas along the 
major axis of the galaxy, associated with enhanced [N II] 
A6584 emission. In addition, this data also provides ev¬ 
idence for a jet-induced outflow with embedded ionisa¬ 
tion cone along the minor axis of the galaxy. The high 
velocity dispersion along the major axis may be due to 
a variety of causes. First, it could be simply line of sight 
overlap between the two expanding bubbles. Second, it 
could be due to a physical interaction between the am¬ 
bient medium and the radio jets. Third, it could be due 
to intrinsic velocity dispersion and shocks in a turbulent 
accretion flow. In all of these scenarios, shocks play a 
major role, and the presence of thermal X-ray emission 
co-extensive with the region of high velocity dispersion 
tends to support this picture. The temperature of this X- 
ray plasma {kT = 0.4— 0.5 keV), if shock-heated, would 
require shocks with a velocity of ^ 400 km s“^. 

Figure reveals the presence of dusty (accretion?) 
streams as well as both clumpy and filamentary emis¬ 
sion line gas in the region occupied by the thermal X-ray 
plasma. Figure I shows that the [N II] 6584A emission 
is particularly strong in this region, which is bounded by 
the X-shaped region of high [O III] 5007A emission which 
coincide with the inner boundary of the expanding large- 
scale bubbles of gas seen in Figured 

We are thus led to a model which nas both inflows and 
outflows on the large scale. An accretion flow is falling 
into a small few arc sec. diameter circum-nuclear region. 
The accretion flow is mostly confined within 45 degrees 
of the major axis of the galaxy. At the same time, a 
hot plasma consisting of both relativistic gas from the 
nucleus and a hot X-ray plasma with kT ~ 0.4 — 0.5 keV 
is rising buoyantly at velocities of at least 150km s“^ 
through the large scale accretion flow. This outflow is 
goverened by the potential of the galaxy such that the 
bubble rises towards the minor axis of the galaxy. The in¬ 
teraction between these rising bubbles of hot plasma and 
the large-scale accretion flow is the fundamental cause of 
the regions of large velocity dispersion seen in the WiFeS 
data. 
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Fig. 6.— The HST FOS spectrum of NGC 1052 (solid black line), the fitted line and continuum spectrum (blue line for the UV or red 
line for optical wavelengths), and the residuals (gray line). For the continuum, a stellar template is used down to 3550A(thin red line), and 
a Chebyshev polynomial fit shortward of that (thin blue line). 


3.2. Accretion Shocks 

In our model, the hot plasma acts as a source of vis¬ 
cosity in the accreting gas, allowing the infalling gas to 
settle at sub-Keplerian velocity. In these circumstances, 
an accretion shock will form on the two faces of the dense 
central accretion disk-like structure. The strongly shock- 
enhanced emission in the accretion shock region pro¬ 
duces the central bright disk so evident in Figure and 
mapped in detail in t he sub- arcse c ond re solution integral 
field observations of Sugai et al. (2005). These authors 
find that this region is characterised by a relatively slow 
rotation with a velocity amplitude ~ 100— I20km s“^ at 


a radial distance of ~ I50pc. The inner outflow axis de¬ 
fined by the radio axis is strongly inclined to the rotation 
axis of this accretion disk by about 50° 

Shocks may also arise internally in the central accretion 
disk if the orbits of the accreting material are somewhat 
ch aotic. Such a inodel w as suggested for the case of M87 
by Dopita et al. (1997). In such a case, radiative inter- 
nal shocks serve to increase the binding energy of the 
accepted material. In other words, these shocks act as 
a source of viscosity in the accreted material, helping to 
convey gas towards the central engine. 

3.3. Cocoon Shocks 


































































Dopita et al. 


The mis-alignment of the radio jet axis and the ac¬ 
cretion disk rotation axis has important consequences 
resulting from the interaction of the radio plasma with 


the (r elatively thick) accretion disk. In the Sugai et al. 

(20051 data, this interaction shows as two [0 ili | AhOO/ 


emission ridges associated with the Kadler et al. (2004a I 
radio jets. The approaching (ENE) jet produces a ridge 
of [O III] emission blue-shifted with respect to the sys¬ 
temic velocity by -400 km s“^ out to a radial distance 
of about 75pc, while the receding radio jet gives a red- 
shifted [O III] emission ridge with a radial velocity of 
200km s“^ out to a radial distance of ^ lOOpc. Beyond 
this radius, the mean radial velocity of both ridges falls 
rapidly towards the systemic velocity. At the same time, 
the [O III]/H/3 ratio rises rapidly. 

Here we interpret these observations in terms of the 
cocoon shock model for the intera ction of a radio jet 


with an ambient interstellar medium (Bicknell et al. 1997 
19981. In this model, the radio jet passes though a ter- 


mmation shock, and becomes surrounded by a cocoon 
of hot plasma which is a mixture of relativistic plasma 
and thermal plasma entrained or swept up by the jet 
material. It is this cocoon which gives rise to diffuse 
radio emission and possibly some of the diffuse X-rays. 
The optical emission arises not from this component, but 
from a strong radiative wall shock propagating into the 
surrounding medium as a result of the overpressure of 
the cocoon material. These radiative shocks may also 
give rise to thermal X-ray emission if the shock velocity 
is high enough (> 200 km s“^). 

In the case of NGC 1052, the mis-alignment of the jet 
and accretion disk axes ensures that these radiative wall 
shocks sit in the already shocked material of the accretion 
disk and are much brighter on the side of the jet closest 
to the mid plane of th e accretion disk. Fro m the observed 
radial velocities from Sugai et al. (20051, we infer that 
the cocoon shock velocity is of order ~ 300 km s“^. If 
the direction of jet propagation is almost in the plane of 
the sky, then the way in which the radial velocity falls to 
the systemic velocity at the head of the radio jets can be 
explained in terms of projection. On the other hand, the 
increase of the [O III]/H/3 ratio can be ascribed to lower 
density, faster, and/or partially-radiative shocks at the 
head of the cocoon. 

The geometry of this model is schematically illustrated 
for clarity in Figure in which we choose a viewing angle 
close to the mid-plane of the galaxy. The mis-alignment 
of the inner accretion disk and the resulting one-sided na¬ 
ture of the cocoon shock is evident, as is the escape of the 
relativistic plasma into the expanding buoyant bubbles 
of gas surrounded by weak optical emission lines seen in 
the WiFeS observations, Figure 

In our cocoon shock model, the radiative wall shock 
is propagating into a medium which has already been 
shocked to high density by the accretion flow. Thus, the 
initial density is orders of magnitude greater than the 
initial density in the accretion shocks. This two-shock 
model then would provide a natural explanation for the 
high- and low-density components observed to be present 
in the FOS observations. 

We now proceed to develop this model in more detail. 
For this purpose we use the photoionisation and shock 
modelling code MAPPINGS v4.2 (Sutherland et ah, in 


preparation). This code now incorporates a global re¬ 
vision of the input atomic data for all elements up to 
and including Zinc, as well as a full revision of the non¬ 
equilibrium ionisation/cooling solver. As such, it is more 
comprehensive from an atomic physics point of view, as 
well as being faster and more stable than previous ver¬ 
sions of the code. 

4. MODELLING THE SPECTRUM OF NGC 1052 
4.1. The Accretion Shock 

All our modelling is done with a dusty plasma having a 
3x solar abundance set and depletion factors as defined 


Dopita et al. (2013). Such high abundances are ap- 


propriate for the centr es of massive and e volved galaxies 


in the local universe (Zahid et al. 2013). We set up a 
medium in which the component ot magnetic pressure 
support and the gas pressure are equal (equipartition). 
The component of the magnetic field transverse to the di¬ 
rection of propagation of the shocks limits the maximum 
compression of the post-shock gas that can be achieved. 
Defining the Alfven Mach Number AIa as the ratio of 
the shock velocity to the Alfven velocity va = /Air pq 
where B is the transverse component of the magnetic 
field and po is the pre-shock density, the maximum com¬ 
pression factor in the shock, pi/po, is given by 


P2 


( 1 ) 


An accretion shock produced in a medium which is ini¬ 
tially in equipartition is already magnetically-dominated 
in terms of its internal pressure support by the time it 
has cooled enough to emit in the [S II] AA6717,31 lines. 
Thus, the ISM pressure which provides the best fit to 
the density-sensitive [S II] AA6717/6731 ratio provides 
a good estimate of the initial pressure in the infalling 
gas for a given accretion shock velocity. On the basis 
of the observed velocity dispersion, we estimate the ac¬ 
cretion shock velocity to be about 150 km s“^, but in 
order to attempt to constrain it further we investigated 
the effect of changing the accretion shock velocity from 
80 to 260 km s“^. With an accretion shock velocity of 
150 km s“^, the parameters of the pre-shock plasma are 
found to be r = lO^K, nn = 24 cm“^, B — 36/rG, 
log P/k = 6.2. Such a pressure is not unreasonable for 
an Elliptical galaxy. For example, i n a detailed study of 
five nearby relaxed giant Ellipticals Werner et al. (2012) 
finds central pressures (in the hot X-ray plasma) averag- 
ing log P/k = 6.5. 

4.2. The EUV diffuse field 

For the atomic phase, we assume this to be pre-ionised 
by the diffuse X-ray spectrum (the molecular gas is too 
dense for this t o be a significant sou rce of pre-ionisation 
in this phase). Kadler et al. (2004b) gives the absorbed 


luminosity of L — (1.4) x 10"^^ erg s ^ and an (intrinsic) 
luminosity in the X-rays of L = (1.7 —2.0) x 10^^ erg s“^. 
In addition, they find a very flat power-law photon index 
r = 0.2 — 0.3. The measured hydrogen absorption col¬ 
umn density of A^h (6 — 7) x 10^^ cm“^ is consistent 
with a nuclear extinction of at least Ay 3mag. For our 
input radiation field we have adopted a theoretical spec¬ 
trum for a black hole of Mbh = IO^Mq accreting at an 
Eddington factor of 10“^ taken from the library of Done 
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Fig. 7.— A schematic picture of the accretion flow, accretion shock and interacting cocoon shock model for NGC 1052 presented in this 

paper (figure courtesy of Dr. David INichollsj. 
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vided by the [O III] lines. As the energy density in 
the radiation field comes to dominate over the energy 
density in the shocks, and the ionisation parameter 
exceeds about logC/ = —3.0, the [O III] /H/3 ratio 
increases to eventually approach the values typical of 
Seyfert galaxies. At the same time, as the ionisation 
parameter increases, the temperature in the [O III] zone 
decreases so that the JO III] AA4363/5007 ratio decreases 
as the shock component of the spectrum decreases in 
relative importance. This is illustrated in Figure for 
the case Vs = 170 km s“^. In this figure, we have used 
the mixing fraction of the coc oon shock defined by the 
observations (see Section |4.4[ below) The best fit to 
observation is obtained with Fx = 0.22 erg cm“^. 

Since the ionisation parameter characterising the X- 
ray photon field is determined by the geometry of the 
diffuse medium with respect to the AGN - which we do 
cannot exactly determine a priori, we treat the dilution 
factor of the X-ray field as a (somewhat) free parameter. 


In the full family of models we have used three values of 
the X-ray field corresponding to Fx = 0.06, 0.22 and 0.66 
erg cm“^, which is achieved at radial distances of 2.8 x 
10^°, 2.2 X 10^° and 1.3 x 10^° cm from the central engine, 
respectively. This should be compared with the size of 
the central bright emission region in the Ha emission, 
which is ro ughly 3 x 10^° cm in radius at the distance of 
NGC 1052 dSugai et al.|[200^ . 


4.3. The Cocoon Shock 

The observations constrain the velocity of the cocoon 
shock to be between 200 and 300 km s“^. Since the co¬ 
coon shock is propagating into the accretion disk material 
which has already been processed through the accretion 
shock, we set the pre-shock pressure equal to post-shock 
pressure in the accretion shock, and once again assume 
equipartition. We also set the ionisation state in equilib¬ 
rium with the EUV radiation field. Due to the low ionisa¬ 
tion parameter of the local radiation field, the pre-shock 
gas is only partly ionised; H+/H = 0.7. This gives a set of 
pre-shock parameters: Te = 3000K, nn = 28000 cm“^, 
B = 600/rG, and logP/fc = 8.2 cm“^K. This very high 
pressure would have to be comparable to the pressure in 
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Fig. 8.— The variation of the ratio of the theoretical to observed 
[O III] /H/3 and [O III] AA4363/5007 ratios as a function of the 
ionising flux in the EUV radiation field. The curves are shown for 
a shock velocity of Vs = 170 km s~^ and a flux contribution from 
the molecular component of 0.36. These curves illustrate how tight 
limits can be placed on the intensity of the radiation field using the 
observed [O III] excitation. 


the radio plasma which is driving the cocoon shock. As a 
result of the high pre-shock density, the electron density 
in the cooling zones of the shock exceeds 10® cm“®, and 
consequently both the cooling time and cooling length 
in the shock are very short, of the order of 2 years and 
3 X 10^"^cm, respectively. 


4.4. Constraining the Cocoon Shock Fraction, cj) 

In our model, lines such as [S II] AA6717, 6731, [O 11] 
AA3727,3729, [N II] AA6548,6584 and to a lesser ex¬ 
tent [O III] AA4959, 5007 have their origin predominantly 
in the accretion shock. On the other hand, lines pro¬ 
duced in dense gas such as [S II] AA4069,4076, [O II] 
AA2470, [N II] AA5755 and [O III] AA4363 as well as 
other lines with a high critical density such as [O I] 
AA6300, 6363 are produced in the cocoon shock compo¬ 
nent. This difference provides the means of determining 
the relative emission fraction, (p, of the cocoon shocks 
in the combined spectrum as seen in the FOS aper¬ 
ture. The total flux is defined in terms of the separate 
fluxes from the cocoon shock and the accretion shock as 

Ftot. ^Acocoon 4“ (1 0 ) ^accretion ■ 

This difference provides a very direct way of constrain¬ 
ing (j). In Figure]^ we show the sensitivity of five ratios 
to the shock velocity in the accretion flow and mixing 
fraction of the cocoon shock component. The mixing 
fraction is very tightly constrained, and found to lie be¬ 
tween 0.28 < < 0.35, regardless of the shock velocity 

in the accretion shock which is here varied from 80 up 
to 200 km s“^. We cannot direct tly estimate the m ixing 
fraction from the observations of Sugai et al. (2005), but 
these are at least consistent with the tact that the cocoon 
shock does not dominate within the FOS aperture. 


4.5. Best Fit Models 

The best fit model is derived from a global fit to the 
emission line spectrum. For this purpose we compare 
the modulus of the logarithmic difference between the¬ 
ory and the observations of NGC 1052 for the predicted 
line intensities (relative to H/3), separately for the UV 
and the optical lines, and attempt to minimise this dif¬ 
ference by varying the shock velocity while keeping the 
molecular fraction, UV radiation field and pre-shock gas 
pressure constant. This procedure is known mathemati¬ 
cally as minimising the LI norm in the logarithmic space 
□ and it allows us to fit both bright lines and faint lines 
simultaneously to similar goodness of fit. 

The result of a run with fixed cocoon shock velocity 
(300 km s“^) and varying accretion shock velocity is 
shown in Figure Marginally we obtain the best-fit 
model for an accretion shock velocity of 150 km s“^. 
However, what is most striking about this figure is the 
insensitivity of the solution to accretion shock velocity - 
essentially no value in the computed range is dehnitely 
excluded. 

The result is also insensitive to the assumed cocoon 
shock velocity. All that we require is that the ram pres¬ 
sure of the shock is sufficient to compress the gas up to 
the range implied by the observed line ratios which are 
sensitive to high density. For illustration of this point, 
we include in Table 1 four models which provide an al¬ 
most equally good fits to the observations. Model A and 
Model B both have an EUV field Fx = 0.07 erg cm“^ 
and an accretion shock velocity of 150 km s“^, but differ 
in their cocoon shock velocities (200 and 300 km s“^, 
respectively). Model C and D have a stronger EUV 
field; 0.22 erg cm“^, but have the same accretion ve¬ 
locity (150 km s“^). They likewise differ in their cocoon 
shock velocities (200 and 300 km s“^, respectively). 

Note that all of these models provide a good descrip¬ 
tion to the full range of ionisation states observed, rang¬ 
ing from O I (IP = 13.61eV) up to Ne V (IP = 126.2eV). 
A very good fit is obtained in the optical, with the ex¬ 
ception of the (N I] line - which we find to be an excep¬ 
tionally difficult line to model correctly as it arises in a 
thin zone where hydrogen transitions from the ionised to 
neutral state. It is worth remarking that the intensity of 
the He II A4686 line predicted by the models and given 
in Table 1 also looks to be in good agreement w i th the 
actual measurement of this line by Sugai et al. (2005) 
(their Figure 1). 


Comparing our results with the Gabel et al. 
2-phase dusty photoionisation model, we hnd t 


( |M| 

rat our 


model gives an appreciably improved fit to the UV lines. 
This is because the high ionisation zones of our models 
are appreciably hotter than in photoionisation models, 
thanks to the influence of shock heating. In general the 
UV lines are enhanced in our model. Likewise the [Ne V] 
lines in the near UV are mu ch stronger in our model. In 
the red, |Gabel et al. (2000) note that their models un- 
derestimate the ]U I AA6300,6363, [N II] AA6548,6584 
and [S II] AA67I7,6731 lines. Doubtless the fit would 
have been improved had they adopted the higher chem¬ 
ical abundances of these elements such as we have used 
in our model. Clearly a reasonable ht to the observed 


S00I http://rorasa.wordpress.com/2012/05/13/10-norm-ll-nonn-12-norra-l-infinity-norm/ 





























Multiple shocks in NGC 1052 


11 


1.00 

1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 

1.00 

1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 
“ [S II] vs. [0 II] 

pr 0.50 

- [S II] vs. [S II] 

.Q 

O 

CD 

^ 0.50 

0.6 

N 

CO 

CO 

N 

CO 

'o) 0.00 

o 

: r c o^cc 

0.0 0.2 0.4 

£0 

O 

CJ5 

° 0.00 

0.4 ^ 

: / ^ : 

^ -0.50 

Mixing Fraction 

S' 

N 

C\J 

1-0.50 

./ 

0.0 

CO 

N 

CO 

CJ) 

o 

-1.00 

1 ■ ■ ■ ■ 1 ' ' ■ ■ 1 ■ ■ ■ ■ 1 ' ' ' ■ 1 

N 

CJ) 

O 

-1.00 

r 

1 ■ ■ ■ ■ 1 ' ' ' ' 1 ■ ■ ■ ■ 1 ■ ' ' ' 1 


-1.0 -0.5 0.0 0.5 1.0 

log(4068,76/6717,31 ),h:log{4068,76/6717,31 )at,s. 


-1.0 -0.5 0.0 0.5 1.0 

log(4068,76/6717,31 ),f,.-log(4068,76/6717,31 )obs. 




Iog(4068,76/6717,31 ),b.-log{4068,76/6717,31 > 055 . Iog(4068,76/6717,31 ),b. -Iog(4068,76/6717,31 )obs. 


Fig. 9.— The variation of the density-sensitive ratios of [S II], [O II], JN II] and [O III] as a function of accretion shock velocity and 
mixing fraction of the cocoon shock component for a cocoon shock velocity of 200 km s“^. We show the ratio of the theoretical line ratio 
to that observed, so by definition, all the observational points are at the origin of the coordinates. The errors for the observational points 
are assumed to given by the photometric measurement errors. The accretion shock velocity is poorly constrained in these figures, since it 
varies over the range 80-200 km s“^ for each curve of fixed mixing fraction, but the mixing fraction itself is closely constrained by these 
data in the range 0.32 it 0.04. 


spectrum can be had with a 2-component photoionisa¬ 
tion model but it would still need to be made consis¬ 
tent with the observed spatial and dynamical structure 
of NGC 1052. 

5. DISCUSSION 
5.1. Black Hole Mass 

The model d eveloped in the pre vious section and the 
observations of Sugai et al. (2005) can provide interest¬ 
ing insights on both the Black hole mass and the power 
in the radio jet of NGC 1052. Let us first consider what 
constraint can be placed on the mass of t he black hole 
in NG C 1052 from the dynamical data. [Sugai et al.| 
(2005) quote the amplitude of the rotation velocity as 
100-130 km s“^ at a radius of 145 pc (1.5 arcsec). The 
inclination of the accretion disk at this radius can only 
be crudely estimated as ~ 30 — 50° from the shape of the 
line-free continuum distribution shown in their Figure 2, 
or the HST Ha observation of Figure Assuming that 
the gravitational potential this close to the nucleus is 


dominated by the black hole, we (very crudely) estimate 
Mbh = (4- 14) X 10®Mq. 

An alternative approach is offered by the measure¬ 
ment of the central velocity dispersion in the WiFeS 
data of Figure 230 km s“^. Using the black hole 
mass: velocity d i spersi on relationship i nitially derived by 
Trema ine et al.j (12002 1 and refined by Graham & Scott 
l 2013| would then imply a black hole mass in the range 
-^BH = (2 —7) X 1O®M0. This is probably be the more re¬ 
liable estimate of the two, although both are very crude. 
Nonetheless, they at least give an idea of the size of the 
central engine. 


5.2. Power of the Radio Jet 

The power of the radio jet can be estimated from th e 
jet - cocoon theory developed by||Bicknell et al.j (|1997|). 
In this, the evolution is driven by the jet parameters, the 
jet energy flux Ej, and the relativistic (3 = vjc. Define 
the r—axis to be in the direction along the jets, and the 
z—axis to be the perpendicular direction. If the mean 
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Fig. 10.— The mean absolute logarithmic difference between the 
line fluxes predicted by theory and the observations of NGC 1052. 
These are here plotted as a function of accretion shock velocity 
for a flxed cocoon shock velocity of 300 km s“^. For the UV we 
have fitted simultaneously to 14 lines, while for the optical the fit 
is for 22 emission lines. The accretion shock velocity is only weakly 
constrained, but for shock velocities below 130 km s“^ the fit is 
appreciably worse in the UV. This figure illustrates the insensitivity 
of our models to the actual accretion shock velocity. Our models 
are also insensitive to varying cocoon shock velocity. 

pressure in the cocoon is P, then towards the head of 
the cocoon (in the r direction) the “head” pressure in the 
direction of propagation is higher than the mean pressure 
in the cocoon. This is caused by the relatively small 
area, A, over which the jet termination shock is jittering 
and the stagnation of the cocoon flow in the direction of 
propagation. Here we can take the pressure as a factor 
C times the average lobe pressure where, typically ( ^ 
2 - 10. However, given that the jet appears to be more 
directly impinging on the cocoon of NGC 1052, we will 
assume C "^1 for the purpose of calculation. The velocity 
of advance of the jet, and the velocity of the wall shocks 
surrounding the cavity excavated by the jet are given by, 
respectively: 


dr 

dt 



1/2 

^ ( 1/2 

P' 

pcA 


[p\ 


and 


dz 

dt 


- 11/2 


1/2 


( 2 ) 


( 3 ) 


where p is the density in the surrounding pre-shock 
medium, and Ej, /3, and A have been defined above. 
The 300 km s“^ cocoon shock model, the post-shock 
pressure, P is found to be 5.2 x 10“® dyne cm“^, and 
the pre-shock density, p = 7.8 x 10“^° gm cm“^. The 
relativistic /3 = 0.38 has been measured from the proper 
motion of hot-spots in the j ets at at a rad i o wav elengths 
in the MOJAVE survey of Lister et al. (2013). From 
these observations, evidence tor precession m the jet is 


inferred. All that remains is to estimate A. Noting that 
the jet is highly collimated we assume a jet opening an¬ 
gle of < 5 degrees, which implies an area A < 10^® cm^. 
Substituting these into equation]^ gives us Ej < 4 x 10"^"^ 
erg s“^. This jet energy flux is comparable to those in- 


ferred for the GPS and CSS radio sources by |Bicknell| 
eFallpWTl). 


5.3. Mass Accretion Rate 

Let us now test whether the accretion shock arises 
solely at an accretion shock into the outer accretion disk 
or has in addition contributions from internal dissipative 
shocks in the accr etion structure itself , as suggested for 
the case of M87 by|Dopita et al.|(|1997D. The key param¬ 
eter which the model should be capable of reproducing is 
the observed absolute H/3 flux, which is a direct measure 
of the rate of increase of binding energy of the accreted 
gas. From the de-reddened flux given in Table 1 and 
the distance of 19.5 Mpc, we estimate Lh /3 4 x 10®® 
erg s“®. The theoretical surface flux in H/3 in our best fit 
model is S'h /3 = (2 —4) x 10 “®erg cn: i~®s~®. The radius of 


the bright accretion structure from Sugai et al. (2005) is 
about 100 pc. Therefore, the total area of the accretion 
disk (top and bottom) is A ^ 5 x 10^® cm®, we predict a 
luminosity of L/j ^ (1 — 2 ) x 10 ®®, in fairly good agree¬ 
ment with what is observed. We conclude that internal 
dissipative shocks within the accretion disk itself cannot 
contribute more than 50% to the observed luminosity. 

In the infall model the accretion rate is M = Apvs, 
where p is the mean density of the infalling gas; ^ 
5 X 10®® gm cm“®, and Vs 150 km s“® is the accretion 
shock velocity. This gives M ~ 6 M 0 yr“®. This accre¬ 
tion rate seems remarkably high. However, we can check 
this by estimating the column density in the accreting 
material, using this to estimate the extinction it would 
produce, and comparing this with what is observed. For 
simplicity we assume that the accretion velocity remains 
constant, so that the accretion flow density falls off as 
r“®. The hydrogen column density is therefore (approx¬ 
imately) the radius of the accretion disk times the pre¬ 
shock density nn ~ 24 cm“®. This implies a column 
density of ~ 7 x 10®® cm“®. If the gas to dust r atio is 
similar to that of the Galaxy (Bohlin et al. 1978), then 
Ay ~ 3, which is a little higher than the Ay ~ (1.1 —1.4) 
inferred from the emis sion line spectrum of the central 
disk (see Section 2.2.2) but is probably consistent with 
the observations when the clumpy nature of the accretion 
flow is taken into account. 

Although the accretion rate is high, only a small frac¬ 
tion of this accreted matter can find its way down to 
feed the central Black Hole. The rad io jet parameters of 
NGC 1052 estimated in Section 5.2 imply (for a Black 
Hole mass of (2 — 7) x IO^Mq) an Eddington fraction 
of < (0.025 — 0.1) and an accretion rate onto the Black 
Hole of < 10“® M 0 yr“®. Thus, either the radio jet is 
very effective in ejecting the accreted gas, or we are see¬ 
ing NGC 1052 in a special epoch in which the accretion 
rate into the outer accretion disk is very high, while the 
accretion from the inner accretion disk to the Black Hole 
is very low. 

5.4. Comparison with other LINERs 
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Although we have shown that the optical emission lu¬ 
minosity in NGC 1052 most likely results from accre¬ 
tion and mechanical luminosity in the jet. However, we 
should ask whether this scenario is valid also for other 
LINERs, particularly those hosted by spiral galaxies, 
where radio luminosities and radio jets are not so com¬ 
mon as in elliptical galaxies, and where the accretion 
rates could be lower. 

We have WiFeS data on two other galaxies, NGC 6812 
and NGC 7213, taken from the S7 first data release (Do- 
pita et al. 2015). These LINER galaxies have many 
features in common with NGC 1052 and cover a range of 
types. While NGC 1052 itself is an E4 galaxy, NGC 3031 
is of SA(s)ab type, NGC 6812 is classified as an SO, and 
NGC 7213 is SA(s)a type in the de Vaucouleurs classifi¬ 
cation scheme. The remarkable similarity of their spectra 
(extracted from the nuclear region with an aperture of 
4 arcsec in diameter) is shown in Figure [TT] NG C 7213 
has been cl assified as a LINER or a Seyfm 1.5 (| Evans 
et al.||1996| based upon the broad Ha visible in Figure 


11 


owever, it is clear that a broad-line LINER classi- 


hcation is more appropriate. 



One of the best-studied LINERs in an spiral galaxy 
host is NGC 3031, better known as M81. This also con¬ 
tains a flat-spectr um core, with clear evid ence for a pro¬ 
cessing radio jet (Martf-Vidal et al. 2011), which would 
set up ideal conditions for the type of jet-accretion disk 
interaction proposed here. M81 co ntains a black hole 


with a mass of about 2 x 1O’’M0 (Martf-Vidal et al. 


2011) and reference s therein), which displays variable 


hard X-ray activity ( ]Pian et al.||2010 ). In the soft X-ray 
component there is clear evidence for thermal plasma 
components with two well-constrained temperatures of 
0.18F0.04 keV and 0.64F0.04 keV ( |Page et al.|20'Mj ). If 
produced by shocks, this would require shock velocities 
of 250 and 470 km s“^, respectively. 

M81, like NGC 1052, contains a compact (~ 140pc di- 
a meter) disk of gas bright in [N II] A 6584 and Ha i mage d 
by |Devereux et al.| (|1997|) - see also |Pogge et al.| (|2000|). 
Leading into this disk is a spiral-like dust lan e, i ' ne U V 
spe ctrum of this galax y has been studied by Ho (19961 
and Maoz et al. ( 1998|). These spectra show many of the 
same features as NGC 1052, but with more evidence of 
an underlying BLR seen in Mg II A2800 and in Ha and 
H/3 and in the strong featureless continuum. We also ob¬ 
tained high signal to noise optical-UV spectroscopic data 
under project GO-6532 (investigators: Dopita & Ko- 
ratkar) using the Faint Object Spectrograph (EOS). This 
data is presented in Figure [T^ in the same format as Fig¬ 
ure]^ The similarity of the two spectra is remarkable. In 
particular, we see the same evidence of a two-component 
density structure. For example, the S II AA6731/6717 
ratio in NGC 3031 is indicative of an electron density 
Tie ~ 1200, while the S II AA4069,76/6717,31 ratio is 


characteristic of densities ~ 3 x 10®. 

Regrettably, the presence of the broad-line component 
makes the determination of the reddening effectively im¬ 
possible in this relatively low-resolution EOS spectrum, 
as we cannot reliably ascribe a Balmer decrement to 
the LINER gas. However, on the assumption that the 
dust extinction is negligible, we can estimate the density- 
sensitive line ratios, given in Table 2. These indicate den¬ 
sities of roughly a factor of two greater than NGC 1052 
in both components of density. In the two-shock model 
advocated here, this is most likely accommodated by an 
increase in the mixing fraction of the cocoon shock (co¬ 
coon shock emission to total emission) to 0.5 — 0.6 - c.f. 
Fig[^ presumably thanks to the lower attenuation of this 
component. Increasing the dust extinction would drive 
this estimate to even higher values. The key point how¬ 
ever is that M81 - like NGC 1052 - requires a two-density 
structure which is naturally explained by the two-shock 
model proposed here. 


Fig. 11. — The WiFeS spectra of three LINER galaxies extracted 
from the nuclear region with an aperture of 4 arcsec. These LIN¬ 
ERs are remarkably similar in their emission line spectrum, despite 
the different host galaxy types. 


All of these galax ies - like NGC 105 2 


are flat spec- 
These exhibit 


trum radio sources (Healey et al. 2007) 
significant polarisation and are beiieved to show appre¬ 
ciable Doppler boosting related to a strong jet compo¬ 
nent. In particular, NGC 7213 conta ins a compact and 
variable radio core (Blank et al.||2005 ). This galaxy, like 
NGC 1052, also possesses an amorpnous core emission 
region ~ 130 pc in diameter. 


6. CONCLUSIONS 

In this paper, we have investigated the complex dy¬ 
namical structure of the prototypical LINER NGC 1052. 
We find that on the large scale (using the WiFeS instru¬ 
ment) (^ 20 arcsec), we have a pair of buoyant bubbles 
of plasma rising into and shocking the extended inter¬ 
stellar medium of the galaxies along the minor axis of 
the galaxy. A portion of these bubbles shows enhanced 
[O III]/H/3 line ratios, and appears to be photo ionised 
from the central engine, while we also see evidence for 
a turbulent accretion flow into a plane roughly aligned 
with the major axis of the galaxy, and in the sense of 

















































14 


Dopita et al. 


■^oo-T- 2-^ CD 

^Sln^^CDCD 

5>>>s>>>>>>>^^ii 

WOOOWOWOW WOOZ lOO 


^ ^ h' O) 
CO O) o o 

a> a> S S 


CO^U^<U ,_^0O ^ CO OJCM iS5 

O)co mmioiolo JG^j-cd h~ oo 

^coco CO CO [it 

T-'T-COCOCOCOCOCO^C\jC\jC\j^^'«!t'«!t 

CVJCVJ CVJCVJOJOJCVJ ^ ^ >>0000^2 = = = = = = 

= == = = = =— <]j<i) <i) <i) a)<u==<]j<]j <]j<i)<i) 

ZZ OOOOOWLi-Li. Li. Li. ZZOOLi.Li. Li.Li.Li. Li.S2 


C\J OJ 




I I 


V 


W/ \l 



Fig. 12.— The HST FOS spectrum of NGC 3031 (solid black line), the fitted line and continuum spectrum (blue line for the UV or 
red line for optical wavelengths), and the residuals (gray line). Compare with Figure^for NGC 1052. The spectra are very similar in all 
salient features, except that we see a broad line spectrum in NGC 3031 seen clearly inrhe Mg II A2800 and the Ha lines. 


rotation of the stars. With HST imaging on scales of 
^2 — 4 arcsec we see a complex dusty accretion flow and 
a bright c entral disk. This dis k is resolved by the obser¬ 
vations of Sugai et al. (20051 into a rotating turbulent 
disk and a fast bipolar outflow which is well-aligned with 
the inner radio jet, but mis-aligned with the ~ 1.0 arcsec 
disk. As a consequence, the jet interacts directly with its 
own accretion disk. 

We have built a self-consistent 2-shock model consist¬ 
ing of accretion shocks from an ISM with log P/k ^ 
6.2K cm“^ located at the surface of the accretion disk, 
and a pair of incomplete cocoon shocks which are being 
driven into the accretion disk by the ram pressure of the 


radio jet. With such a model, we can explain the ob¬ 
served velocity dispersion, and the full UV-optical emis¬ 
sion line spectrum of the nuclear region of NGC 1052 
obtained with FOS on the HST. The two shock model 
naturally explains why there seem to be two components 
of very different density, one with rie 10^ and the other 
with Ue > 10® cm“®. We find that the theoretical spec¬ 
trum is very sensitive to the ratio of the cocoon shock 
emission to the accretion shock emission, but that it is 
very insensitive to the assumed shock velocities, either 
in the accretion flow or in the jet-shocked cocoon. Thus, 
this type of a model may well be applicable generically to 
other members of the LINER family of AGN. If so, much 
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of the optical emission luminosity of this class would be 
the result of accretion and mechanical luminosity in the 
jet, rather than due to EUV radiation from the central 
black hole. 
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TABLE 1 

Extracted nuclear line fluxes for NGC 1052 (de-reddened by Ay = 1.05) compared 
WITH the models DESCRIBED IN SECTION 4.5. 


A (A) 


Line ID 

Obs. Fluxi Model A^ 

(H/5 = 100) 

Model B® 

Model 

Model D® 

1548,51® 


CIV 

72 ± I 2 FI 

174 

119 

144 

123 

1640 


He II 

74 ± 9 '-' 

72 

69 

73 

79 

1661,66 


[O III] 

30 ± 10 

67 

57 

70 

60 

1883,92 


Si III] 

51 ± 7 

12 

10 

13 

12 

1907,09 


C III] 

129 ±7 

69 

67 

82 

92 

2139,43 


N II] 

38 ±5 

8 

8 

8 

8 

2323,35 


C II], Si II] 

107± 13 

116 

134 

116 

122 

2470 


[O II] 

41 ± 3 

58 

54 

82 

89 

2422,24 


JNe IV] 

37 ± 4 

30 

27 

28 

24 

2785,93 


Mg II 

67 ± 6 

55 

66 

45 

48 

3345 


JNe V] 

4± 2 

5 

5 

5 

4 

3426 


Ne V 

12 ±3 

13 

12 

15 

12 

3726,29 


[O II] 

403 ±8 

595 

617 

665 

647 

3869 


]Ne III] 

76 ±8 

59 

65 

71 

97 

3889 


He I,H I 

26 ±2 

26 

28 

26 

27 

3968,70 


]Ne III], He 

39 ±8 

34 

36 

35 

43 

4069,76 


[S II] 

77 ±5 

109 

121 

84 

80 

4100 


H<5 

19 ±3 

26 

26 

26 

26 

4340 


H7 

50 ±3 

46 

47 

47 

47 

4363 


JO III] 

21 ± 3 

16 

13 

16 

15 

4686 


He II 


10 

10 

10 

11 

4861 


H/3 

100 ±3 

100 

100 

100 

100 

4959 


[O III] 

85 ±3 

73 

64 

78 

76 

5007 


[O III] 

241 ±4 

211 

186 

225 

220 

5198,200 


[N I] 

31 ± 4 

2 

3 

2 

3 

5755 


JN II] 

12 ±5 

13 

14 

15 

17 

5876 


He I] 

17 ± 6 

20 

22 

20 

21 

6300 


[O I] 

183 ±5 

183 

241 

121 

175 

6363 


[O I] 

54 ±3 

58 

77 

39 

56 

6548 


[N II] 

111 ± 7 

110 

110 

126 

123 

6563 


Ha 

274 ± 20 

294 

294 

294 

294 

6584 


[N II] 

330 ± 21 

322 

324 

372 

362 

6716 


[S II] 

143 ±5 

190 

188 

184 

186 

6731 


[S II] 

131 ± 5 

169 

178 

195 

190 

^ Absolute de-reddened flux in aperture: F{Hq) - 

= 9.04 X 10“'^"' (erg cm~ 



2 Fx = 0.07 

erg 

cm~^, accretion shock velocity 150 km s“^, 

cocoon shock 

velocity 200 km s“^. 

® Fx = 0.07 

erg 

cm“^, accretion shock velocity 150 km s“^, 

cocoon shock 

velocity 300 km s“^. 

* F^ = 0.22 

erg 

cm“^, accretion shock velocity 150 km s~^, 

cocoon shock 

velocity 200 km s~^. 

® Fx = 0.22 

erg 

cm~^, accretion shock velocity 150 km s~^, 

cocoon shock 

velocity 300 km s~^. 

° Line doublets 

are indicated thus (i.e) 1548,51 = 

A1548 -h A1551 



^Errors given 

are measurement errors only. 






Errors due to reddening correction are not computed 


TABLE 2 

A COMPARISON OF THE DENSITY-SENSITIVE LINE RATIOS MEASURED WITH FOS IN NGC 1052 AND IN NGC 3031. 


Line Ratio Ion NGC 1052 NGC 3031 


6731/6717 

ns 

nr 

0.91 

± 

0.06 

1.32 

± 

0.06 

4069,76/6717,31 

fs 

II] 

0.28 

± 

0.05 

0.53 

± 

0.09 

4363/5007 

fo 

III] 

0.087 

± 

0.012 

0.149 

± 

0.049 

2470/3726,9 

fo 

II] 

0.107 

± 

0.009 

0.30 

± 

0.05 

5755/6584 

[N 

II] 

0.027 

± 

0.013 

0.107 

± 

0.011 









